Citations:
Land use/land cover data: http://geoportal.hawaii.gov/datasets/land-use-land-cover-lulc Watershed data: http://geoportal.hawaii.gov/datasets/watersheds
### Wrangling and Exploration:
# Look at LULC:
plot(lu_lc)
#unique(lu_lc$landcover)
# 28 different land use types. One of them is just "0" -- does that mean no data available?
# Could group land use types into categories: Agriculture, Residential, Urban, Aquatic?
#### TYPES ####
#[1] "Cropland and Pasture"
#[2] "Commercial and Services"
#[3] "Residential"
#[4] "Evergreen Forest Land"
#[5] "Other Urban or Built-up Land"
#[6] "Mixed Rangeland"
#[7] "Industrial"
#[8] "Streams and Canals"
#[9] "Orchards, Groves, Vineyards, Nurseries and Ornamental Horticultural Areas"
#[10] "Shrub and Brush Rangeland"
#[11] "Forested Wetland"
#[12] "Reservoirs"
#[13] "Nonforested Wetland"
#[14] "Bare Exposed Rock"
#[15] "Sandy Areas Other than Beaches"
#[16] "Transportation, Communications and Utilities"
#[17] "Herbaceous Rangeland"
#[18] "Beaches"
#[19] "Other Agricultural Land"
#[20] "Lakes"
#[21] "Strip Mines, Quarries, and Gravel Pits"
#[22] "Mixed Barren Land"
#[23] "Bays and Estuaries"
#[24] "Mixed Urban or Built-up Land"
#[25] "Transitional Areas"
#[26] "0"
#[27] "Industrial and Commercial Complexes"
#[28] "Confined Feeding Operations"
lu_lc_categories <- lu_lc %>%
mutate(
landcover = case_when (
landcover %in% c("Industrial and Commercial Complexes",
"Mixed Urban or Built-up Land",
"Transportation, Communications and Utilities",
"Industrial",
"Commercial and Services",
"Other Urban or Built-up Land" ) ~ "Urban", #Group this into a broader "Urban" category
landcover %in% c("Confined Feeding Operations",
"Other Agricultural Land",
"Orchards, Groves, Vineyards, Nurseries and Ornamental Horticultural Areas",
"Cropland and Pasture") ~ "Agricultural", # Group these into a broader "Agricultral" category
landcover %in% c("Bays and Estuaries",
"Lakes",
"Reservoirs",
"Streams and Canals") ~ "Waterways", # Group these into a broader "Waterways" category
landcover %in% c("Nonforested Wetland",
"Forested Wetland") ~ "Wetlands", # "Wetlands" category
landcover %in% c("Mixed Rangeland",
"Shrub and Brush Rangeland",
"Herbaceous Rangeland") ~ "Rangeland",
TRUE ~ landcover)) %>% # Keep everything else as is
filter(landcover != 0) # Filter out the "0" landcover (no data?)
unique(lu_lc_categories$landcover)
## [1] "Agricultural"
## [2] "Urban"
## [3] "Residential"
## [4] "Evergreen Forest Land"
## [5] "Rangeland"
## [6] "Waterways"
## [7] "Wetlands"
## [8] "Bare Exposed Rock"
## [9] "Sandy Areas Other than Beaches"
## [10] "Beaches"
## [11] "Strip Mines, Quarries, and Gravel Pits"
## [12] "Mixed Barren Land"
## [13] "Transitional Areas"
# Now 13 categories, little bit better. Could group more if necessary?
plot(lu_lc_categories)
# Still have area and perimeter. Not necessary but could be useful in visualizing the total area of each land use type.
### Maybe only visualize certain land uses??
### Visualize Land Use / Land Cover for just Big Island
lulc_plot <- ggplot(data = lu_lc_categories) +
geom_sf(data = lu_lc_categories,
color = NA, # Gets rid of gray coloring borders which made categories difficult to see
aes(fill = landcover)) +
scale_fill_paletteer_d("ggsci::springfield_simpsons") +
xlab("Longitude") +
ylab("Latitude") +
ggtitle("Land Use and Land Cover on Big Island, Hawaii") +
labs(fill = "Landcover") +
theme_minimal() +
coord_sf(xlim = c(-156.5, -154.5), ylim = c(18.5, 20.5), expand = FALSE) + # Clip to just Big Island
theme(panel.background = element_rect(fill = "paleturquoise")) # Add background light blue color to make it look like ocean. Is it too bright? Need another color maybe.
## Could add more details in this plot? Maybe add ocean color background. Could also zoom in to one island? Is it too difficult to see individual land uses? Or maybe need to decrease categories.
lulc_plot
Figure 1. Land Use and Land Cover on Big Island, Hawaii. Land use and land cover for thirteen different categories.
Make it interactive!
## Make an interactive map using Tmap
tmap_mode("view") # Tmap set to interactive viewing
lulc_tmap <- tm_basemap("Esri.WorldImagery") + # Add general basemap
tm_shape(lu_lc_categories) +
tm_fill(col = "landcover", palette = "Set1") + # Fill by landuse
tm_legend(legend.position = c("right", "bottom")) +
tm_layout(title = "Land Use and Land Cover in Hawaii")
#leaflet::providers$Stamen.Watercolor
lulc_tmap